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1 Introduction 

We continue our test of the Master Constraint Programme [T| for Loop Quantum Gravity (LQG) 
[HJ III IE] which we started in the companion paper 2 and will continue in jHHUEl- The Master 
Constraint Programme is a new idea to improve on the current situation with the Hamiltonian 
constraint operator for LQG ||]. In short, progress on the solution of the Hamiltonian constraint 
has been slow because of a technical reason: the Hamiltonian constraints themselves are not 
spatially diffeomorphism invariant. This means that one cannot first solve the spatial diffeo- 
morphism constraints and then the Hamiltonian constraints because the latter do not preserve 
the space of solutions to the spatial diffeomorphism constraint On the other hand, the 

space of solutions to the spatial diffeomorphism constraint |1U| is relatively easy to construct 
starting from the spatially diffeomorphism invariant representations on which LQG is based 
which are therefore very natural to use and, moreover, essentially unique. Therefore one would 
really like to keep these structures. The Master Constraint Programme removes that techni- 
cal obstacle by replacing the Hamiltonian constraints by a single Master Constraint which is a 
spatially diffeomorphism invariant integral of squares of the individual Hamiltonian constraints 
which encodes all the necessary information about the constraint surface and the associated 
invariants. See e.g. |2j for a full discussion of these issues. Notice that the idea of squaring 
constraints is not new, see e.g. |12) . however, our concrete implementation is new and also the 
Direct Integral Decomposition (DID) method for solving them, see p] |2j for all the details. 

The Master Constraint for four dimensional General Relativity will appear in ^1] but before 
we test its semiclassical limit, e.g. using the methods of ^3 and try to solve it by DID 
methods we want to test the programme in the series of papers (21 El IH El • We begin by studying 
finite dimensional systems, in particular: a finite number of Abelean constraints linear in the 
momenta which will also play an important role for a system with second class constraints, 
first class constraints with structure constants at most quadratic in the momenta and first class 
constraints linear momenta with structure functions. 
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2 Finite Number of Abelean First Class Constraints Linear in 
the Momenta 



Our first example is a system with configuration manifold K™, coordinatized by x l , i = 1, . . . , n 
and m < n commuting constraints 

Ci=Pi i = l,...,m , (2.1) 

where the p^s are the conjugated momenta to the x l, s. 

All phase space functions which do not depend on Xi are Dirac observables, i.e. functions 
which commute with the constraints (on the constraint hypersurface). A Dirac observable which 
depends on pi, i = 1, . . . , m is equivalent to the Dirac observable obtained from the first one by 
setting pi = (since these two observables will coincide on the constraint hypersurface). There- 
fore it is sufficient to consider observables which are independent of the first m configuration 
observables and of the first m conjugated momenta. A canonical choise for an observable algebra 
is the one generated by x l ,pi, i = (m + 1), . . . , n. 

To quantize the system, we will start with an auxilary Hilbert space ^(R™) on which the 
operators x l act as multiplication operators and the momenta as derivatives, i.e. we use the 
standard Schrodinger representation: 

^V(x)=x^(x) Pi^(x) = -ihdi^(x) . (2.2) 

According to the Master Constraint Programme, we have to consider the spectral resolution 
of the Master Constraint 

m 

M = -h 2 Y,df=:-h 2 A m . (2.3) 
i=i 

This spectral resolution can be constructed with the help of the spectral resolutions of the 
operators pi = —ihdi, which are well known to be given by the Fourier transform: 

V(x)=/ <6 fei ,V(x 1 ,...,x i - 1 ,-,x i+1 ,x") >ib ki (x l )dk t (2.4) 

with generalized eigenfunctions 

< b ki ,tp(x 1 ,...,x i - 1 ,-,x i+1 ,x n ) >i= / ■^=exp(-ik i x i )ij(x)dx i (2.5) 

J IR. 

and ip G -Z^R"). The spectrum of pi is therefore spec(j5j) = R. Moreover the spectral projectors 
of pi and pj commute, so we can achieve a simultaneous diagonalization of all the C,'s (that is 
the Fourier transform with respect to (x 1 , . . . , x m ): 



with 



x)= / <nM-)>^--- 

[Y[b ki (x i ))dk 1 ---dk m (2.6) 



< Ylhi(-)^{--- ,x m+1 ,...,x n ) > {m) = / ( \[b ki {x l )\il>{x)dx l ■ ■ ■ dx rn . (2.7) 
i=i jRm ^i=i ' 

This decomposition of a function ip G L2(R n ) corresponds to a decomposition of the Hilbert 
space L2(R n ) into a direct integral of Hilbert spaces 

L 2 (R n )^ [ H {kl _ km) dk 1 ...dk m . (2.8) 
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where each T~t(ki,...,k m ) 1S isomorphic to L 2 (]R™~ m ). 

Since the Master Constraint Operator is a polynomial of the Cj's, according to spectral 
calculus we already achieved the spectral resolution of the Master Constraint. The spectrum of 
M is given by spec(M) = clos{^™ 1 kf, fcj G 1} = JL|_. The generalized eigenfunctions of M to 
the (generalized) eigenvalue can be read off from (|2.6f) to be 



1 m 

*(x) = ^=^exp(^^)| fci=0 ^(x m+1 ,-.- ,x n ) 

\ ' i=l 

-^{x m+1 ,--- ,x n ) with ^eL 2 (Sl n - m ) (2.£ 



(V2tt) 

i.e. functions which do not depend on (x , . . . ,x m ). So one could conclude that the physical 
Hilbert space is Hm q) — L2(R n ~ m ), a space of functions of (x m+1 , . . . ,x n ). The action of the 
elementary Dirac observables x l and pi for i = (m + 1), . . . ,n is well defined on this physical 
Hilbert space. 

However we would like to go through the explicit procedure for constructing the direct 
integral decomposition for a seperable Hilbert space, since we will use it later on in another 
example. The general procedure is explained in detail in |2j. 

We will work in the fourier-transformed picture, i.e. the Hilbert space is L^iJ^ 1 ) as a function 
space of the momenta (pi, . . . ,p n ) and the Master Constraint Operator becomes the multiplica- 
tion operator M = Yli^xPi- 

To begin with we have to choose a set of orthonormal vectors from which H = L2(R n ) 

can be generated through repeated applications of the Master Constraint Operator M. To this 
end we use the fact (see (U), that 

m m 

{exp ( - ± ^ pj) ( Pj) k k e N, I e Q (2.10) 

3=1 3=1 

is a basis for L2(R m ). Here ip s i are harmonic polynomials 1 of degree s on ]H m , and the index 
I takes values in some finite index set l s , which depends on m and s. Obviously the set (|2.1U|) 
can be generated by applying repeatedly M (more correctly M where M = M <g> ^L 2 {'K m ~ n )^ so 
that M is an operator on L2(R m ')) from the set 

m 

{^:=exp(-i^p2)V s ,/|sGN,/G[ s } . (2.11) 

3=1 

Moreover the spaces generated from different f2 Si ; are mutually orthogonal. (This can be seen 
if one introduces spherical coordinates. Then the restrictions of the harmonic polynomials to 
the unit sphere, which are given by the spherical harmonics, are orthogonal in L,2(S m ~ 1 , duj n ^\) 
where duj n -\ is the uniform measure on the sphere S m ~ 1 .) So we conclude that the set of 
orthonormal vectors 

m 

{fi sZt :=iV sZ exp(-i^p|)V^(8)^(Pm+l,---,Pn)|sGN,Ze[ S! tef} . (2.12) 

3=1 

where N s i are normalization constants and {(frt}tet is an orthonormal basis of L2(R n_m ) fulfills 
all the demands required above, i.e. 



L 2 (R n ) = ^2 (B span{p(M)n slt \p polynomial} . (2.13) 

s,l,t 



1 i.e. homogeneous polynomials annihilated by the Laplace operator on W 
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Now the vectors £l s n are tensor products of the form $l s [ t = £l s i ® cj> t where Q s i £ H\ := L2(]R m ) 
and (pt £ Ti-2 '■= L2(R n ~ m ) and furthermore M = M ® acts only on the first factor. Hence 
we have for the spectral measures fi s it(\): 

Hsu -.= <n s i®(f> t ,(@(\- M)n s i)®(f> t > H =<n s i,e(\- M)n s i> Hl \\(f> t \\ H2 

= <n sl ,e(\-M)n sl > ni =:fi sl , (2.14) 

so that we only need to consider the spectral measures fj, s i with respect to the operator M 
on 7ii. Additionally, using the rotational symmetry of the Master Constraint Operator (and 
an idea outlined in |2J|) one can simplify the calculations even more: The space of harmonic 
polynomials on Jt m of degree s is an irreducible module for the rotation group 0(m) under the 
left regular representation, which acts as 

U(R) : ^j{p) i — > ^{R- 1 ■ p) (2.15) 

on the space of functions over R m . Here R is a rotation matrix for Il m . This representation is 
unitary if considered on the Hilbert space Ti\ = ^(JL" 1 ). The irreducibilty of the representation 
ensures that one can generate all vectors {Q s i}i&i by applying rotations U(R) to just one vector, 
say f2 s := ^I s q. But these rotations leave the spectral projectors of M invariant (since M 
commutes with the rotations U (R) ) , so that we have 

fi sl = <n sl ,e(\-M)n sl >=<u(R)n s0 ,e(x-M)u(R)n s0 > 

= <n s0 ,e(\- M)tt s0 >= n s0 =: n s (2.16) 

for an appropriate rotation R. Here one can see, that it is very helpful to know the symmetries 
of the Master Constraint Operator, i.e. unitary operators commuting with M. These may come 
from (exponentiated) strong Dirac observables, that is operators which commute with all the 
constraints on the whole Hilbert space H. Examples for these are operators of the type 1-^ 1 <g> U2, 
i.e. unitary operators which act only on TC2 = L2(Il"~ m ). We used this kind of symmetry in 
Q2.14|) . However the C/(i?)'s from above are not of this type: The classical counterparts of their 
generators (i.e. angular momentum on lH m ) Poisson-commutes with the Master Constraint but 
they also vanish on the constraint hypersurface {pi = 0, % < m}. Nevertheless the rotation group 
0(m) is useful in constructing the direct Hilbert space decomposition. Later we will see, that 
the vanishing of its (classical) generators on the constraint hypersurface corresponds to the fact, 
that the representation of the rotation group on the induced Hilbert space is trivial. 

So we just need to select for each s € N a particular homogeneous harmonic polynomial ip s 
of degree s. We choose 

ip s (pi,---,Pm) = (Pi+ip2) s ; s = 0,l,2,... (2.17) 
so that the vectors O s become 



n s = = _ p-2W|p^ (2.18) 



1 _i.a 

where we introduced the coordinates r m = YliLi Pi > r i = an d the angle <~p = arctan(pi/j>2 
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With this at hand we can compute the spectral measures [a s (we will assume that m > 4): 

m-(a) = < n s ,@(\-M)n s > Hl 

1 



' m — '2~' m-2 



7T m ' Z Si Jjcn 

m ) 2 ; / / dr 2 / do; m -3 / dr m _ 2 9(A - r^r^+V; 
7r m / 2 s! Jo Jk + J5™-3 Jr + 



>+l_m-3 -rj, 
'm-2 e 



r 2 =r m cos<£ 7T m / 2 S! 

»"m-2=»* m Sill^ 



T #/ dr m G(A-r^+ m - 1 (cos ( / ) ) 2s+1 (sin0r-V 



2nVo\(S m ~ 3 ) s! r(m/2 - 1) 1 
y =rl W 2 5 ! 2r(s + m/2) 2 7 R+ ' 

tv j 7^ / ^e(A- y ) y s+ W2-i e ^ (219) 
r(s + m/2) y R+ 

From the second to the third line we transformed the cartesian coordinates (p3,---,Pm) to 
spherical coordinates on R m ~ 2 with radial coordinate r 2 n _ 2 = YUL^Pi an( ^ an gels varying over 
5<m-3. Vol(5 m ~ 3 ) is the euclidian volume of the (m — 3)-dimensional unit sphere. From the 
third to the fourth line we again introduced spherical coordinates (r m ,<f>) for the two radii r 2 
and r m _2, so that = r\ + r 2 n _ 2 . Since r 2 ,r m _ 2 are positive, (j) varies over [0, vr/2]. Then we 
integrated over eft and used a new integration variable y = for the remaining integral. In the 
last step we used that Vol(S m - 3 ) = 2^ m / 2 -7r(m/2 - 1). 

To obtain the final spectral measure we have to sum over all the measures fi s n = fi s 
multiplied by constants a s i t - We choose the constants such that Y^it a sit = 2 _s_1 , so that the 
sum over all the a s ;t's is one. The spectral measure becomes 

oo 

m(a) = a ^ vsuw = 2 ~ s ~ v.(a) 

s,l,t s=0 

_ ayi ^/ *„_ 2 e ( A-^)rl.«c_-| e -5. 

s=0 ''^ 

_ 7TVol(S™- 3 ) f ,f ^ _qa J* ir^ m -3„-^ 



7T 



m/2 



/ dr 2 / dr m _ 2 e(A-r^)r 2 e2 r 2 r ™_ 3 e- 



7rV0l( ^" ) / cfa/ cir m „2e(A-r 2 a _2-2x)e-e-^CZ 3 



z=r 2 2 /2 W 2 



7rVol(S m - 3 ) f ^ oM _ 2 ^ _I^I„2 , 

7T m / 2 



dr m _ 2 9(A - ri_ 2 )(l - e -2 A +2^- 2 ) e -^-2r r m - 3 



m-2 

1 /",„,. 



y=rl_ 2 r(m/2-l) 



/ dy&(X-y)(l - e -5 x+ h) e -yy m / 2 - 2 (2.20) 



Here we inserted for /i s the third line of (|2.19|) . exchanged integration and summation, performed 
a variable transformation x = r 2 /2 in the fourth line, and integrated over the new variable x in 
the fifth line. Finally we changed to the integration variable y = r 2 n _ 2 in the last line and used 
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the explicit expression 27r m / 2 1 /T(m/2 — 1) for Vo\(S m 3 ). This gives for the derivative of //(A): 

d\ T(m/2 - 1) y R+ 

+ 1 / di,Pi(\ - 7,Uip-5 A +^^-f 1 , m /2-2 

r(m/2 - 1) A + 



^3T)^ dye(A-y)(i e -^) e -V 

2rW2-i) e ^ A l^ e(A ^ )e "^ m/2 " 2 • (2 ' 21) 



Next we have to calculate the Radon-Nikodym derivatives p s = p s i t of p s with respect to p. 
Since both measures are absolutely continuous with respect to the Lebesgue measure, we can 
write 



= dps(X) = dp s (X)/d\ = 2r(m/2- l) Afl+m/2 _ lc _| A r ' 



dp(X) djjL(\)/d\ r(m/2 + s) 



dye 2 y y 

o 



y m/2-2 



(2.22) 



For < A < oo the derivatives p s (X) will be some strictly positive numbers for all s 6 N. But 
for the limit A^Owe apply L'hospital's rule, resulting in 

2r(m/2-l) (s + m/2 - l)A s+m /2~2 _ i^+m/2-i 
lim p s (A) = — ^ ; f lim — ? 2.23 

A^o P v ; r(m/2 + s) A^o A" 1 / 2 - 2 V 7 

which gives p s (0) = for s > and po(0) = 2. 

Now we can construct the induced Hilbert space 7Y®(0). We notize that there is only one 
linearly independent harmonic polynomial of degree zero, namely V'oo = 1- Therefore 7Y®(0) has 
an orthonormal basis {e t }t<=t which corresponds to the set of vectors {^ootjtet- Hence we can 
identify H®(0) with the Hilbert space L2(R"~ m ) of functions in the variables (p m +i, . . . ,p n )- 
Interestingly, because of p s it{X) > for A > 0, the induced Hilbert spaces Ti®[X) for A > are 
in some sense much bigger, since they have a basis {e s ;t | s 6 N, I G I, t G t} corresponding to 
all the vectors Q s i t . The latter can be seen as a basis in L2(S m ~ 1 , \~^~ duj m -i) (g) ^(R™ - ™), 
therefore we can identify 7i®{\) with L 2 (S m - 1 , A^^-i) ® L 2 (R n ^ m ) for A > 0. 

We would like to mention, that the space of (formal) solutions to the Master Constraint 
Operator is much bigger than the space of functions, which are independent of the first m 
coordinates. For instance, functions, harmonic in the first m coordinates, are solutions of the 
master constraint but they are unphysical because they do not solve the individual constraints Cj. 
These solutions automatically do not appear in the spectral resolution of the Master Constraint. 
The intuitive reason for this will be discussed in the conclusions. 



3 A Second Class System 

Here we will discuss a simple second class system, given by the constraints 

Cj = x l and Bi = pi i = 1, . . . , m 

(3.1) 

on the phase space R 2ra . (The x l denote configuration variables and the pi their conjugated 
momenta.) The Poisson brackets are given by {Ci,Bj} = 5ij and all other Poisson brackets 
vanish. Because of the Heisenberg Uncertainty relations one cannot expect to find eigenfunctions 
of the Master Constraint Operator corresponding to the eigenvalue zero. We will therefore alter 
the Master Constraint and validate whether this gives sensible results. 

A complete set of observables is given by phase space functions which are independent on 
the first m configuration variables and momenta. 
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As in the last case we will quantize the system by choosing the auxilary Hilbert space £ 2 (R n ) 
of functions of the configuration variables ?/>(x). The operators x 1 act as multiplication operators 
and the momentum operators Pi as derivatives. 

The Master Constraint Operator is defined as the sum of the squares of the constraints but 
we can as well consider the following slight variation of this prescription: 

m 

M=^^+ U , 2 (ff) , (3.2) 

i=i 

where uj{ are m positive constants. This Master Constraint coincides with the Hamiltonian for 
m uncoupled harmonic oscillators with different frequencies. This Hamiltonian has pure point 
and positive spectrum, the lowest eigenvalue being \ = | Ya=i 

Since zero is not in the spectrum of the Master Constraint Operator we have to alter M, such 
that its spectrum includes zero. Of course we have to check in the end, whether this procedure 
gives a sensible quantum theory. The simplest thing one can do, is to subtract Ao from the 
Master Constraint to obtain M = M — Ao- This is equivalent to the normal ordering of the 
Master Constraint Operator. This could be done also for the limit m — > oo if we multiply the 
individual operators in the sum of 1)3. 2 jl by positive constants Q m with ^2 m Q m < oo. This 
condition will naturally reappear in the free field theory examples of [I]. 

Now the spectrum of M includes zero and we can construct the induced Hilbert space to the 
eigenvalue zero. To this end we have to find a cyclic basis, which we will choose to be a tensor 
basis of the following kind: 

^ni,...,n m ,k = /rai,...,ra m (^1) • • • j ^m) ® fofc^m+l , • • • , X n ) (^-3) 

where / ni ,...,n m > n % £ N are the (normalized) eigenfunctions of M (considered on L2(R m )) and 
{hk | k £ W} is an orthonormal basis in L2(R n_m ). The application of M to the vectors Q ni , ...,n m ,k 
just multiplies them with the corresponding eigenvalue so that we need them all in our cyclic 
basis. The associated spectral measures are 

A*ni,...,n ra ,fc(A) =< ttn u ...,n m ,k, ©(A - M )^ ni ,...,n m ,k >= ©(A - A niv .., nm ) (3.4) 

where A niv .. inm = | YslLi n i^i- Only measures with m = = ■ ■ ■ = n m = have the point 
zero in their support. The final spectral measure can be defined to be 

A* (A) = Yl 2 ~ {ni+1) ■ ■ ■ 2" (nm+1) • 2-( fc+1 Vn 1 ,...,n m , fc (A) = ©(A) + „(A) (3.5) 

ni,...,n m ,fcGlN 

where u{\) = Y, ni ,...,n m >i,km 2" (ni+1) • • • 2-^+^ ■2-^ k+1 ^ nit ...^ k {\) is a pure point measure 
which does not have support at zero. The relevant Radon-Nikodym derivatives for the induced 
Hilbert space to the eigenvalue zero are therefore 

Pni=o,...,n ro =o,fc(0) = 1 , (3.6) 

hence we can interpret the set {^m=o,...,n m =o,fc I k 6 N} as a basis in the induced Hilbert space. 
So as expected for a pure point spectrum the induced Hilbert space is just the (proper) eigenspace 
to the eigenvalue zero. This eigenspace can be identified with the space L2(R ra_m ) by mapping 
^ni=o,...,n m =o,A; to the basis vector hf c (x m+ i, . . . ,x n ) in L2(R n ~ m ). Also the action of Dirac 
observables (that is quantized phase space functions, which do not depend on the first m config- 
uration variables and momenta) coincides on the null eigenspace of M and on L2(R n ~ m ). We 
can therefore conclude, that the physical Hilbert space is given by L2(R n_m ), which is the result 
one would expect beforehand. 



8 



4 Finite Number of First Class Constraints at most Quadratic 
in the Momenta 



If the constraints generate a (semi-simple) compact Lie group it is in general straightforward to 
apply the Master Constraint Programme. The Master Constraint Operator coincides in this case 
with the Casimir of the Lie group and has a pure point spectrum. The direct integral decompo- 
sition of the kinematical Hilbert space (which in this case is truly a direct sum decomposition) 
is equivalent to the reduction of the given representation of the Lie group into irreducible rep- 
resentations and the physical Hilbert space corresponds to the isotypical component 2 of the 
equivalence class of the trivial representation. 

4.1 SU(2) Model with Compact Gauge Orbits 

Here we consider the configuration space R 3 with the three so(3)-generators as constraints: 

U = e%x j p k , {L u Lj} = e%L k (4.1) 

where e\j = is totally antisymmetric with €123 = 1 and we summed over repeated in- 
dices. (In the following, indices will be raised and lowered with respect to the metric git = 
diag(+l,+l,+l).) 

The observable algebra of this system is generated by 

d = x l pi e + = x l Xi e~ = p l pi 

{d,e ± } = ^2e ± {e + ,e~} = Ad . (4.2) 

It constitutes an sl(2, H)-algebra. We have the identity 

-d 2 + e + e" = LiV (4.3) 

between the Casimirs of the constraint and observable algebra. 

4.1.1 Quantization 

We start with the auxilary Hilbert space £ 2 (R 3 ) of square integrable functions of the coordinates. 
The momentum operators are pj = —i{h)dj and the x 3 act as multiplication operators. There 
arises no factor ordering ambiguity for the quantization of the constraints, but for the observable 
algebra to close, we have to choose: 

d = \{x l pi +f>iX l ) = x % pi - \%h e + = x l xi e~ = fipi (4.4) 

The commutators between constraints and between observables are then obtained by replacing 
the Poisson bracket with •] . 
The identity (|4.3j) is altered to: 

d 2 - i(e+r + e~e + ) - \h 2 = -UV . (4.5) 

For the implementation of the Master Constraint Progamme we have to construct the direct 
integral decomposition with respect to the Master Constraint Operator 

M := L\ + l\ + L\ = - (d 2 - l(e + e~ + e"e+) - f h 2 ^j (4.6) 

2 Compact Lie groups are completely reducible. The isotypical component of a given equivalence class within 
a reducible representation is the direct sum of all its irreducible representations into which it can be decomposed 
and which lie in the given equialence class. 
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The Master Constraint Operator is the Casimir of SO(3) on £ 2 (R 3 ). Its spectrum and its 
(normalized) eigenfunction are well known, the latter are given by the spherical harmonics. To 
discuss these we make a coordinate transformation to spherical coordinates: 

xi = rcos(^)sin(0) (j) G [0, 2ir) and 9 G [0, n) 
X2 = r sin (</>) sin (0) 
X3 = rcos(9) 

(4.7) 

d 3 x = r 2 sin(9)dr d9 d(fr 

3 = Ll + Ll + Ll = -h 2 \^-dl + -±-d e (smedej\ (4.8) 

Lsm 6 v smfl J 

p 2 = -n 2 A = -ft 2 (^d r (r 2 d r ) + ^d e (sin9d s ) + -^d; 

x 2 = r 2 (4.9) 

The eigenfunctions of M = L 2 to the eigenvalue h 2 1(1+1), I G N are of the form Yi m (9, </>) R(r) 
where Yj m , I G N, — / < m < I are the spherical harmonics on the two-dimensional sphere S 2 
and R(r) is an arbitrary function in L2(R+, r 2 dr). To discuss the direct integral decomposition 
of the kinematical Hilbert space L2(Il 3 ) we have to find a cyclic basis of this Hilbert space. 
We will choose the eigenbasis of H := e~ + e + , i.e. the three-dimensional harmonic oscillator 
Hamiltonian. Its (normalized) eigenfunctions are given by |22j 

2 

Vtam(r, e, 4>) = N nl r l exp(-^r) M(-n, I + |, \r 2 ) Y lm {9, 0) (4.10) 

where the index n takes values in N, the constant N n i is a normalization constant and M (— n, 1 + 

3 

2' 

-n'of n ro nmnnQV "l-no f 1 1 n rdinn c A/7 / ti / _l 

2' 



- r 2 ) are confluent hyper geometric functions. Since the first argument of these is a whole 
negative number, the functions M(—n, r 2 ) reduce to polynomials in r 2 of order n. Therefore, 
for fixed indices m, I the functions ipnim are polynomials in r with minimal degree r . (From this, 
one can see, that L2(R 3 ) is not equivalent to a tensor product L2CR+, r 2 dr) (g) L2(S 2 , sin 9d9d<p). 
If this would be the case the {i/j n im \ n £ N} should span the whole Hilbert space L2(IH+, r 2 dr) for 
arbitrary fixed indices m, I. However, they span only the subspace of polynomials with minimal 
degree I. We will come back to this point later on.) 

Applying M to the functions ip n i m just multiplies them with h 2 l(l + 1), hence we need them 
all as cyclic vectors Cl n i m . The associated spectral measures can easily be calculated to be 

^nim(A) =< Ipnim, 9(A - M^nlm > = < Ipnlm, ©(A - H 2 l(l + l))lpnlm >= ©(A - h 2 l(l + 1)) . 

(4.11) 

The spectral measures are pure point, they are homogeneous in n and m and have only the 
point ft? 1(1 + 1) as support. Hence the kinematical Hilbert space decomposes into eigenspaces 
of M in the following way 

00 

L 2 (R 3 ) = ^2 clos(span{^ nim | n G N, -I < m < I}) . (4.12) 
1=0 

The induced Hilbert space to the eigenvalue A = is given by clos(span{-i/> n oo | n G N}). These 
functions are constant in 9, <j) and are polynomials in r 2 of order n (weighted with a gaussian 
factor). The above set can be taken as an orthonormal basis in L2(R+, r 2 dr), therefore we can 
identify the physical Hilbert space with L2(Jl+, r 2 dr). 
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On this (physical) Hilbert space the observable algebra is given by 



e+ = r 2 e~ = -h 2 ^d r (r 2 d r ) d = -ih(rd r + %) (4.13) 



and because of (|4,5|) we have 



d 2 - i(e + e~ + e~e + ) = \h 2 . (4.14) 



4.2 Model with Structure Functions rather than Structure Constants Linear 
in the Momenta 

As an example for a model with structure functions we will discuss a sort of a deformed SO(3) 
constraint algebra 

Ci = x 2 pf - x%p 2 C 2 = x%p 1 -x 1 pf C 3 = x lP2 -x 2 pi (4.15) 

{C 1 ,C 2 } = mnp™- 1 x^- 1 C 3 {C 2 ,C 3 } = C 1 {C 3 ,d} = C 2 (4.16) 

where n and m are positive natural numbers. For n = m = 1 we recover the S'0(3)-algebra. 
One can do a similar deformation for SO(p, g)-algebras. 

The Dirac observables (see |18| ) for this system generate an sZ(2,R) algebra: 

+ 2,2, ^ l-m n+1 - 2,2, * 1— n m+1 

e = x x +x 2 H p H Xo e =Pi+p 2 J * Pq 

1 2 m+n A 6 1 2 m+n s 6 

2 

^ = xi P i+x 2 p 2 -\ ■ — x 3 p 3 (4.17) 

m + n 

{d ) e ± } = +2e ± {e + ,e"} = 4d . (4.18) 

For general m, n these observables commute only on the constraint hypersurface with the con- 
straints. 

We define the Master Constraint as 

M = C\ + C\ + Ci . (4.19) 

The constraints C\ and C 2 do not strongly Poisson-commute with M (for n / 1 or ra / 1), but 
C 3 does. For the rest of this chapter, we will consider the case m = 1 and n odd. For these 
parameter values the phase space function e + is non-negative and commutes on the whole phase 
space R 3 x R 3 with the Master Constraint. This will be helpful for the spectral analysis of the 
quantized Master Constraint. The case n = 1 and m odd can be treated in the same way, using 
Fourier transformation. For m > 1 the observable e + contains negative powers of p 3 , therefore 
the analysis gets more complicated. 

4.2.1 Quantization 

We will quantize this system in the usual way by assigning multiplication operators Xi to the con- 
figuration variables x\ and differential operators pj = —ihdj to the momenta pj. The kinematical 
Hilbert space, we are starting with, is L 2 (R 5 ). 

There arises no factor ordering ambiguity for the quantization of the constraints, and since 
the structure function / = n(x3)™ _1 commutes with C 3 (see (j4.16}0 . it is possible to quantize all 
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the constraints as symmetric operators as we will explain in the conclusion section and to have 
C3 standing to the right of the structure function: 



[C u C 2 ] = ihn{x 3 ) n C 3 [C 2 , C 3 ] = ihCt [C 3 , Ci] = iHC 2 
We have to analyze the Master Constraint Operator 

M = C\ + C\ + C\ , 



(4.20) 



(4.21) 



which is a second order partial differential operator. It can be densely defined and is symmetric 
on the linear span of the Hermite functions. 

To solve the eigenvalue equation for the Master Constraint Operator we will introduce new 
coordinates (t, 9, <p) analogous to the coordinates (r, 9, 4>) in the SO (3) case. However, for n ^ 1 
we have to consider the regions X3 > and x 3 < seperately: For 13 > we define 



and for x 3 < 



xi = tcos(psm8 t>0 0g[O,2tt) 8 £ [0, § ] 

X2 = t sin (j> sin 9 

n + 1 1 2 
x 3 = ( — — (t cos6>)" +1 



^3 = -( ! ^)^(-i COS^^TT G [| )7r ] 



(4.22) 



(4.23) 



n + 3 



The measure is transformed to dx\dx2dx 3 = sin 9 (^j) ™ +1 i ™ +1 1 cos 0| "+ 1 dtdOdtp. With these 
coordinates we have 

(4.24) 



i 2 = e^ 



n+l" 3 

and M does not include derivatives with respect to i. Indeed, for cos 9 > 0, 



M 



1 1 \ zn 

2 1 n -\- L \ n +\ 



{t COS 9) "+ 1 2n-2 



(n — 1 



(icos#) "+ 1 d 2 + ( ^777^ (cos0)«+i 1 r I siu ( 



i sin 



4n 

(t cos 6')™+ 1 
t 2 sin 2 



+ 



n+1 
2 

n+1 



+ 



2n 
n+1 



(4.25) 



and M for cos 9 < is obtained from the above formula by replacing t with —t. The operator 
M simplifies considerably, if we introduce the coordinate 



u = (cos 9) "+ 1 for cos 9 > 



cos6»)™+i for cos#<0 u€[-l,l] . 



(4.26) 



The coordinate u is proportional to the coordinate x 3 and therefore one gets the same operator 
for x 3 positive and for x 3 negative: 



M 



2n-2 

n+1 



(4.27) 



11 1 n+d 

The measure is now dx\dx 2 dx 3 = (' n ^-) n + 1 t n + 1 dtdudcf). 

Since the operators C 3 = —ihd^ and M commute, we can diagonalize them simultaneously. 
Choosing periodic boundary conditions in <p one obtains spec(Cs) = {hk\k G Z} and the null 
eigenspace consists of constant functions in <p. We are interested in the spectrum of M near zero, 
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so it suffices to consider M restricted to the null eigenspace of C3, as it has elsewhere spectrum 
bounded from below by f? (since M = C\ + C 2 + C*|). The restriction of M to this space is 

2 2n-2 

M^ = -n\^)^ 2 ^d u {{l-u^)d u ) , (4.28) 
and can be seen as a product of two commuting operators B\ and B 2 : 

„ 2n-2 

B ^-= h \ r ^)^ 2 i^ ' B 2 :=-dMl-u^)dA . (4.29) 

In the following, we will consider the operator B2 on the Hilbert space L 2 ((—l, l),du) and 
show that its spectrum is purely discrete. Afterwards we will come back to the product of the 
operators B\ ■ B 2 . (The product B\ ■ B 2 cannot be seen as a tensor product of B\ and B 2 since 

q n-f-3 

L 2 (R , dx A ) is not equivalent to L 2 (R + , t n + l dt) (g> L 2 ((— 1, 1), du) tg> L 2 ((0, 2tt), dcfi). Moreover 
one has to be careful, since it is not guaranteed that B 2 has a self- adjoint extension as an 
operator on L 2 (^). This might be the case, because B 2 does not preserve the space of Hermite 
functions in the three variables (x\, x 2 , X3), only the combination B\-B 2 does preserve this space. 
Nevertheless one can consider the operator B 2 defined on the Hilbert space L 2 ((—l, l),du).) 
The operator B 2 is symmetric and positive on the dense domain 

®(B 2 ) = {/ G L 2 ((-l, 1), du)\f G AC(-1, 1), (1 - u n+l )f G AC(-1, 1), 

B 2 -f€ L 2 ((-l, 1), du), lim (1 - u n+l )f{u) = 0} (4.30) 

u—*±l 

where AC denotes the class of absolutely continuous functions. The positivity can be seen by 
using integration by parts 

<f,B 2 -f> L2(u) = j dj{uj{l-u n+1 )d u f{u)du >0 (4.31) 

for functions / G D(B 2 ). A positive and symmetric operator has always selfadjoint extensions 
(see PU). 

The operator B 2 is an ordinary differential operator, more specifically it is a Sturm-Liouville 
operator on the interval I := (—1, 1), so we can utilize the theory of Sturm-Liouville operators, 
see for example |2U| . Sturm-Liouville operators can be classified according to the behaviour of 
eigensolutions at the endpoints of the interval /: A Sturm-Liouville operator A is limit circle at 
the endpoint a if for one A G C all solutions to (A — A)/ = are square integrable near a. One 
can prove, that if this is the case for one A, then it holds for all A G C. It therefore suffices to 
consider the solutions to B 2 ■ f = Xf for one A G C, in particular A = 0. 

One solution for A = is obviously given by fm(u) = 1, the other solution is given by 

Both solutions are square integrable on (—1, 1). The function fy 2 is integrable because for n > 1 
and n odd, the function \ fo 2 (n > 1)| is less than \ fo2( n = 1)1 on the interval (—1,1). But 
fo2( n = 1) = artanh(n) is square integrable on (—1,1). We are therefore in the limit circle case 
for both boundaries u± = ±1. 

If both boundary points are limit circle for an operator A, the following holds (see |2L)|): 
The spectrum of any self adjoint extension is purely discrete, the eigenfunctions are simple 
(multiplicity one) and form an orthonormal basis and the resolvent of A is a Hilbert - Schmidt 
operator. Notice that Hilbert - Schmidt operators are in particular compact, hence the spectrum 
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of the resolvent R Z {A) for z E R not in the spectrum of A has an accumulation point at most at 
zero. Thus A does not have any accumulation point and since M is unbounded we know that 
the eigenvalues actually diverge (when ordered according to size). 

Hence B2 has a discrete simple spectrum, and since foi(u) = 1 is a square integrable solution 
for A = (which fulfills contrary to /02 the boundary conditions), zero is included in the 
spectrum. This holds for any selfadjoint extension of B<i- Now B2 with the domain Q4.3UJ) is 
positive. We choose in the following any selfadjoint extension such that it is still positive and 
such that B\ ■ B2 as an operator in ^(R 3 ) corresponds to the positive self-adjoint extension of 
M, we have chosen before (e.g. the Friedrich extension) 3 . The eigenvalues of B2 are therefore 
positive (that is Xj > and Xj = if and oly if j = 0) . 

Since the operator B2 has a null eigenspace consisting of functions constant in u and therefore 
9, one would expect that the physical Hilbert space consists of functions functionally independent 
of 9 (and eft as was argued above). However we are dealing not just with the operator B2 but 
with the product B\ ■ B2. The operator B\ is a multiplication operator and has continuous 
spectrum on the whole positive axis including zero. 

Therefore we have to discuss the decomposition of the Hilbert space into a direct integral 
Hilbert space and to calculate the spectral measures. To this end we will consider the restriction 
of the Master Constraint Operator to the subspace C3 = as it has elsewhere a spectrum, which 
is bounded from below by ti 2 . Since this subspace coincides with functions in ^(Jt 3 ), which are 
constant in the angle variable (j), we can identify this subspace with 7i' := L2(R+ x ~R,pdpdx3) 
where the variable p is defined by p 2 = x\ + x\ . 

A dense system of vectors in this Hilbert space is generated by the set 

{f pq := p 2 ^x p 3 eM~b 2 ) I W^} • (4-33) 

We need just the even polynomials in p since the variable p has range only in the positive half 
axis 4 R+. Now one has to find a cyclic system of vectors with respect to M. We will begin 
with functions that are polynomials in t 2 = (p 2 + ^iX^ +1 ) weighted with exp(— t 2 ). All these 

functions are annihilated by M, hence to get a cyclic system, we need all the even powers of 
t. (Again, we just need the even powers, since t extends over the positive half axis.) If one 
performs the Gram-Schmidt procedure for this system {i 2fc exp(— t 2 ) \ k G N} one will get an 
ortho-normal set of vectors {vo k | k £ N} which can be identified as a basis in the Hilbert space 

re + 3 

L2(H+ , t "+ 1 dt). The ascociated spectral measures are given by 

p ok (X) =<v 0k ,Q(X- Bf B 2 )v ok >=<v k,&(X- B 1 -0)v ok >=e(X) , (4.34) 

so that the set of vectors {fo/c} is associated with a pure point spectral measure. We will call the 
subspace spanned by these vectors H pp (anticipating that all vectors orthogonal to this subspace 
are associated to spectral measures, which are absolutely continuous). 

Now we have to find the orthogonal complement to TC pp . If one rewrites the functions f pq in 
terms of the coordinates (t, u) one gets 

fpq = p 2q xl = C(p,q)(l - u n+1 ) q u p t 2<1+ ^e-h 2 (4.35) 
where C(p, q) is ap, q dependent constant. It is important to note here that a power of u a is alway 

2a _ 2 

accompanied by a power of at least t n +* . The underlying reason for this is that u ~ x$t "+ 1 . 
Using this fact, it is straightforward to see, that also 

9pq = u p t^ +2q e~^ t2 p, q 6 N (4.36) 

3 The operator Bi is self-adjoint with the domain Dsa(Bi) = {/ £ L 2 {K i )\B 1 ■ f G L 2 (K, 3 )}. 

4 Indeed, polynomials of any positive power of t r could be used in order to construct a basis of the Hilbert space 
by the Gram - Schmidt procedure because any square integrable function t i— * f(t) can be written as t 1— * f r (t r ) 
where fr(s) := /(s 1//r ). Notice that we do not transform the measure to the variable t T here. 
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or 

hpg = P p (u)t^ +2q e-^ t2 (4.37) 

generate a dense subspace in 7i'. Here P P {u) are the Legendre polynomials of order p, that 
is P p {u) is a polynomial with degree p. To see this, notice that (j4.35|) is a polynomial in 
terms of t 2 and x := ut 2 ^ n+1 ' > which explains the statement for g pq . Next notice that P p is a 
polynomial of order p in u which can be written as a polynomial in x of the same order and 
£ 2 (p— *0/( n + 1 ) j = 0, ..,£>. Now any positive power of t can be expanded in terms of the vq^ and 
thus in terms of the t 2q again. 

Since Po(u) is the constant function, the space spanned by {ho q \ q £ N} coincides with the 
space spanned by {vok \k S N}. Moreover P p ,p > is orthogonal to Pq (in L2((— 1, 1), du)), 
so that the orthogonal complement to TL PP is spanned by {h pg |N 9 p > 0,g G I}. We will 
now show, that the spectral measures associated to vectors from this subspace are absolutely 
continuous. In order to do this we have to expand the functions P p (u) £ L2((—l, l),du),p > 
into eigenfunctions of B>2'. 

oo 

P p {u) = Y, a kMu) (4-38) 

fe=i 

where ipk is the normalized k — th eigenfunction of B>2 (and we assume that the corresponding 
eigenvalues A& are ordered and that tpo is the constant function). The constant function (i.e. ipo) 
does not appear in this decomposition since P p is orthogonal (in L2{(— 1, 1), du)) to the constant 
function for p > 0. 

Using this decomposition and abbreviating c n = h 2 (^^-) n + 1 we can write 
Vp'q'pqW ■= <h p > q >,6{\-M)h pq > 

'72. + 1 \ TTXT 



/ / dtdueM-t 2 )t^ +2{ci+q,)+ ^{ YaU k ,(u)) 

/ oo \ 

2n-2 / i — -v „ 



, • I") 

\ fc=i 



T V/ / ^exp(-t 2 )t^ +2 ^') + Sf 
L _i JK+ ^[-1,1] 

2n-2 / 

0(A - c n t ™+i A fc ) a£,a£ ^ fe , (u)V>& («) 
V / dt exp(-t 2 )T 



1 x ' dtexp(-t 2 )t "+i +2(9+?,+ »+i 

fc=l- /R + 

2n — 2 _/ _ 

0(A-c n i— A fc ) a£ a£ (4.39) 



where in the last line we used the ortho- normality of ip k hi -^2(( — 1, 1), du). 

The sum in the last line converges absolutely, since the integral over t can be bounded (k- 
independently) from above by ignoring the ^-function. Then we are left with the sum over the 
absolute values of a p k a?, which can be estimated using the Cauchy-Schwarz inequality and the 
fact that the aP k are the expansion coefficients of P k £ L2((— 1, l),du). 

This shows that all the measures /i P Yp<? are absolutely continuous with respect to the 
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Lebesgue measure on R + . Indeed we have 



- ( n + 1 )^YJ a P( n + l W ( A A ) 3W+ ^ ,)( " +1) 

V 2 J Z.^ a ^ 2CnAfe (n-l)^ eX ^ ( c n A fc j CnAfc^ 

3+p+p , + ( g +g , )(u+l) 

=: A s=i /pY P9 (A) (4.40) 

2n — 2 

where in the second line we performed a coordinate transformation s k = c n Xkt "+ 1 in order to 
solve the delta-function. The sum in the second line converges absolutely (uniformly in A) since 
the exponential factor can be estimated by 1 and the A^ have to be bigger than 1 for some finite 
k because the spectrum of B2 does not have an accumation point as we showed above. Hence 
fp'q'pq is wen defined, in particular f p ' q ' pq (0) > for p' = p (since all terms in the sum are then 
positive). 

To summarize, our Hilbert space 7i' can be decomposed into two subspaces TL PP and TC ac 
where M restricted to H pp has a pure point spectral measure and M restricted to H ac has an 
absolutely continuous spectral measure. We therefore have to discuss the direct integral de- 
composition of 7i pp and 7i ac separately. The 'direct integral decomposition' for Ji pp is straight- 
forward: Since TC PP is a proper eigenspace to to the null eigenvalue of M, the physical Hilbert 
space corresponding to the pure point spectrum coincides with 7i pp (which can be identified 

with L 2 (R+,t^dt)). 

To perform the direct integral Hilbert space decomposition of 7i ac spanned by {h pg |N B p > 
0,g 6 N} we choose an ortho-normal cyclic system {Q m |N 3 m > 0} in that space, such that 
fli = N±hio, where N\ is a normalization constant. Since {h pq |ME5p>0,(/eM} generates H ac 
we can always find coefficients A^ q such that 

00 

«m= £ A p q h Pq • (4-41) 

p=l,9=0 

Consider the spectral measures [i m ,m > 1 associated to this cyclic system 

Mm(A) := <n m ,9(\-M)n m > 

00 

= <n m , 9(X-M) £ A mh Pq > 

p=l,q=0 

00 

= <n m , e(x-M) £ a p>pi> ( 4 - 42 ) 

p=l,q=0 

(P,?)^(1,0) 

since < Q m ,9(X — M)Q\ >= N\ < Q m ,6(X — M)h\o >= due to the defining property of a 
cyclic system. Hence 

00 00 

/UA) = £ £ Vi' A m< h p , q ,,9{X-M)h pq > 

p'=l,q'=Q p=l,9=0 
(P,9)#(l,0) 
00 00 

= E E V q ' A m • ( 4 - 43 ) 

p' = l,q' = p=l,9=0 

(p,g)^(l,0) 
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Differentiating under the sums 5 



dX 2^ 2^ a p'i' pi dX 

p'=l,q'=0 p=l,q=0 
(P,5)^(l,0) 

oo oo . . 

^ ^ -rm Am ■ p+p'+(i+i')(n+i)-3 

p' = l,g'=0 P=l,9=0 

(P,9)^(l,0) 

=: A^<? m (A) (4.44) 

where g m (X) is a non-singular function of A near zero (that is g m (0) is either zero or some finite 
value). This holds for m > 1; for m = 1 we get 



dpi (A) 
dX 



5 



= A^(iV 1 ) 2 /ioio(A) (4.45) 



where /ioio(0) = c> 0. 

The total spectral measure p oc (A) on H ac is given by (if the maximal multiplicity is less than 
infinity, choose appropriate normalization constants different from 2 _m ) 



dy ac (X) l rijm(A) \- d^ m (A) 

dA ' 2 <tt ^, dA 

m=2 

= lA^r(iV 1 ) 2 /ioio(A) + A^T ^2- m 5m (A) . (4.46) 

m=2 

Now we can consider the Radon-Nikodym derivatives in the limit A — > 0. Assume that m > 1: 

6 

t m r dfi rn (X)/dX A«-i# m (A) 

hm p m (A) = hm = hm = = = (4.47) 

dfiacW/dX a^o i A^(iVi)Vioio(A) + A— £ m=2 2- m 5m (A) 

Hence all the p m ,m > 1 vanish at A = 0. But the same calculation for m = gives pi(0) = 2. 
Therefore the induced Hilbert space for A = from H ac is one-dimensional, that is unitar- 
ily equivalent to (D. Putting the contributions from TC PP and H ac together we get TCphys — 

re + 3 

L2(H+, i n+1 ) © <C (where the inner product in (D can be rescaled arbitrarily). The first term cor- 
responds to the (proper) null subspace with respect to M the second to the continious spectrum 
of the multiplication operator in t. 



Remark: 

By the standard theory for Sturm - Liouville operators such as Bi we know that the eigen- 
function ip^ has k zeroes and that the eigenvalues asymptote to A& oc k 2 . By the Weierstrass 
theorem, ip^ can be approximated in the sup - norm arbitrarily well on the compact set [—1,1] 
by polynomials of degree at least k (in order to have k real roots). Thus, the tp^ are actually not 
too different from the standard Legendre polynomials and if they would be really polynomials 
of degree k we could just use the functions Q k = ■p k t 2k ^ n+v > as a cyclic system. Then the above 
calculations would become entirely trivial because the Radon - Nikodym derivative at A = 
would be obviously non - vanishing for the lowest order k = 1. Unfortunately the eigenvectors 
ipk are not polynomials unless n = 1 and so we had to go through this very elaborous analysis. 

5 Strictly speaking one must verify that we may interchange summation and differentiation. This could be done 
for instance by verifying that the series for /x m (A) converges absolutely and uniformly in A at least close to zero 
(the absolute and uniform convergence of the series for jj,(X) then follows). We have not checked that this is the 
case but given the fact that the fj, p i q i pq converge rapidly to zero in the vicinity of A = this should be true. With 
the tools given here, one could check this explicitly by a tedious but straightforward calculation. 
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5 Pure Point and Absolutely Continuous Spectra 



Here we will discuss an example, where, similarly to the previous one, the Master Constraint 
Operator has pure point and absolutely continuous spectrum at zero. The example is sim- 
pler without structure constants and serves the purpose to illustrate an important point when 
choosing the cyclic system. 

We will start with a kinematical Hilbert space L2(R 2 ) and a Master Constraint Operator 

M = (xl + xl)(x 1 p 2 -x 2 pi) 2 (5.1) 

where the xi are multiplication operators and the pi := —ihdi act by differentiation. Apriori one 
would expect that the physical Hilbert space includes the space of functions with zero angular 
momentum and a one-dimensional part which corresponds to the generalized eigenfunction 5{x\)- 
S(x 2 ). 

We will now construct the direct integral decomposition of the kinematical Hilbert space 
with respect to M. First of all, we have to find a cyclic basis for L 2 (R 2 ). To this end we note, 
that L 2 (R 2 ) is spanned by polynomials in x\,x 2 weighted by a gaussian factor. These can be 
generated by the set 

{x^x™ 2 exp(-±(x 2 + x\) | m,n 2 G N} (5.2) 

and also by the set 

{vmn = x™x™L exp(— ix + x_) | n, m G N} (5.3) 
where we defined x± := xi ± ix 2 . In spherical coordinates defined by 

x\ = r cos^ r G R+, <p G [0, 2tt) 

x 2 = r sin0 (5.4) 

the vectors v nm can be expressed as 

Vnm = r n+m exp(^(m - n)) e~ r ^ 2 . (5.5) 

By introducing new indices N := (n — m) G Z and k := \{m + n — \n — m\) G N the set 115. 3 jl 
can also be written as 

{ VNk = r \N\+2k exp ( iiV 0) exp(ir 2 ) | N G Z, k G N} . (5.6) 

Applying repeatedly M to a vector > N ^ generates all the other vectors v^k,^ £ N with 
the same index N. We therefore need just the vno as cyclic vectors. But for N = applying M 
to vok gives zero, so that we also need the whole set {v^k \ k G N}. However, notice that this set 
is not orthonormal, one has therefore to perform a Gram-Schmidt procedure. This will result 
in an orthonormal set of the form {Qok '■= fki 7 " 2 ) ex p( — r2 /2)} where the are polynomials of 
order k in r 2 . This set can be taken as a basis for L 2 (R+,rdr). 

Remark: 

At this point it is important to draw attention to the following subtlety: Recall the definition of 
the tensor product of two Hilbert spaces Hf This is the Hilbert space 7i\®7i 2 consisting of pairs 
ipl (dtp 2 := (ifti,ip 2 ) with ipj G TCj equipped with the inner product < ipi ®ip 2 ,ip'i ®^' 2 >"Hi®w 2 := 
Ili=i ^ ^ji^'j >W 3 - O ne can show that if a£ is a basis for Tij then a£ (g> bff is a basis 
for Hi ®H, 2 . If (Xj,Bj, fJLj) are measure spaces and Tij := L 2 (Xj , dfij) are separable Hilbert 
spaces then one can show by using Fubinis theorem that 7i\ <8> Ti 2 is isometrically isomorphic to 
L 2 (M\ x M 2 , d/ii (8) dfi 2 ) via tp\ ® ip 2 1— > tpi(xi)ip 2 (x 2 ) where ji\ ® \i 2 is the measure on M\ x M 2 
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based on the smallest a— algebra B\ % B2 containing the "rectangles" B\ x _B 2 where Bj E Bj 
and by definition [i\ % ^(-Si x B2) := fJ-i(B '1)1x2^2)- 

Now consider the problem at hand: In polar coordinates we have /^(K 2 , d 2 x) = L2(R 2 , rdrd<p) 
withR 2 = [(R+-{0})x5 1 ]U{(0,0)}. One could now think that since the one point sets {(0, 0)} 
and {0} respectively have d 2 x and dx Lebesgue measure zero respectively that 

L 2 (R 2 ,d 2 x) = L 2 (R 2 - {(0,0)},d 2 x) = L 2 ({R+- {0}) x S 1 ,d 2 x) 

= L 2 ((R + -{0},rdr)®L 2 (S 1 ,d(l>) = L2((Il + ,rdr)®L2(S 1 ,d<t)) (5.7) 

If (|5.7[) would be true then, given some ONB b(l) k (r) for L 2 (R + ,rdr) consisting of polynomials 
in r times e~ r2 ^ 2 obtained via the Gram - Schmidt procedure and an ONB bn\<p) = e in ^/V27r 
for L 2 (S l ,d4>) we would obtain a basis by <g)&„ for L 2 (R 2 ,d 2 x). However, in the tensor product 
we then obtain a dense set of vectors of the form r k e m <t> e - r2 / 2 where the pair (k, n) E No X Z is 
unrestricted. On the other hand the Hermite polynomial basis for L 2 (R 2 ,d 2 x) provide a dense 
set consisting of vectors of the form r \ n \+ 21 e mt t> e r / 2 as we have just seen. We conclude that the 
basis b^ <g> bf^ for L 2 ((R + , rdr) (g) L 2 {S 1 ,d<j)) is overcomplete for L2(R 2 ,d 2 x) and hence these 
Hilbert spaces are not isometrically isomorphic. In other words, the functions r k e m ^> e ~ r I 2 such 
that k — \n\ is not a non - negative and even integer could be expanded in terms of Hermite 
polynomials because they are obviously square integrable. 

What has gone wrong in (|5.7|) is that on IR 2 the coordinates (r, (p) are singular at r = 0. 
The coordinates xi,X2 are globally defined which results in the fact that there is the restriction 
k — \n\ = 21, I = 0,1,2,.. in order that the functions r k e tn ^ are are regular at r = when 
expressed in terms of x\,X2- This topological subtlety that IR 2 is a plane while IH+ x S* 1 is a 
half - infinite cylinder with different orthonormal bases is of outmost importance for the Direct 
Integral Decomposition (DID) because neglecting this difference would result in a much bigger 
physical Hilbert space as we will see shortly. 

The normalization of the vectors vmo gives 

n m = (tt |iV!|)- 1/2 rl iV l exp{iN(f>) exp(-r 2 /2) . (5.8) 
Their associated spectral measure are 

O /'OO 

AtJV0 =<O JV0 ,e(A-M)O A r > = Jj^J o r 2 \ N \exp(-r 2 )e(X-r 2 N 2 )rdr 

- mi ^-n-w-^)* ^ 

which shows, that these measures are absolutely continuous with respect to the Lebesgue mea- 
sure on JL|_. 

The spectral measures associated to the Qok can be calculated to 

Mo* =< V ok , 0(A - M)fl ok >= 9(A) < fl ok , fl ok >= O(A) (5.10) 

which shows, that these measures are of pure point type. 

Hence the kinematical Hilbert space L,2(R 2 ) decomposes into a direct sum of two Hilbert 
spaces TL ac and Tt pp , on which M has either absolutely continuous or pure point spectrum. 7i ac 

^ k 

is defined to be the closure of the span of {M Qnq \ N ^ 0, N G Z, k £ N} and Tt pp is the closure 
of the span of {Oofc I k E N}. We will now discuss the direct integral decomposition on each of 
these two spaces seperately. 

The direct integral decomposition of 7i pp (which can be identified with L2(R+, rdr) is easy 
to obtain, since all the measures no k are the same and have just the point zero in their support. 
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Hence TC pp is alredy decomposed with respect to M - the null eigenspace of M on TC pp coincides 
with 7i pp . The contribution from the pure point part of the spectrum to the physical Hilbert 
space is therefore given by 7i pp . 

For the direct integral decomposition of 7i ac we have to calculate first the total spectral 
measure 



^ac(A) = \ ^2- JV ( mo (A) + /i-7vo(A)) 

JV=1 

roo ( 00 ( r /<2)N \ 

= I I ^^^e(A-iV 2 x)jexp(-x) ( ix . (5.11) 



\N=1 

This gives for the Radon-Nikodym derivative with respect to the Lebesgue measure 

d^ac sr^ 1 / A 



El4 exp(-A/^) . (5.12) 



d\ ^ iV 2 AH \ 2iV 2 

7V=1 v 

The Radon-Nikodym derivative of the measures /i^q with respect to the Lebesgue measure is 

dfim A'-^l 2 

-^r = A7w exp( 7 } ( } 

so that the Radon-Nikodym derivatives of fiNO with respect to [i ac can be calculated to 
pm(X) = p2l£ = xm (iV 2 |iV|))-'exp(-A/^) 

In the limit A — > most of the pnq vanish: 

[2 for IJVI = 1, 
taW= n , , ' 5.15 
A^o I for I TV I > 1. 

Hence the contribution from the absolutely continuous part of the spectrum to the physical 
Hilbert space consists of just two vectors. This contribution corresponds to the fact that the 
delta function 5(xi)5(x2) is a generalized eigenvector of M. 

The total physical Hilbert space is the sum of the contributions from TC pp and TC ac , i-e the 
sum of Tl pp which can be identified with L 2 (R+,rdr) and two vectors {eio,e_io}. Notice that 
if we had made the wrong identification of L2(R 2 , d 2 x) with L2(R+, rdr) ® L>2(S l , d(f>) discussed 
above then we would have had to use the vectors Q' N0 oc e m ^e~ r , n / to get a cyclic 
system. The spectral measures of these vectors all coincide and therefore the corresponding 
Radon - Nikodym derivatives would all be non vanishing at A = and the contribution to the 
physcal Hilbert space from the continuous spectrum would be infinite dimensional which would 
be physically wrong because the constraint x\ + x\ = corresponds to only one point in phase 
space and the correponding physical Hilbert space should be finite dimensional. 



6 Conclusions 

The main lesson learnt in this article is that the Master Constraint Programme can deal essen- 
tially with all situations that one usually encounters when quantizing finite dimensional systems. 
It is even possible to deal with second class constraints as already emphasized by Klauder in 
his Affine Quantization Programme for gravity [T2|. A situation that we have not dealt with 
are constraints which are quadratic in the momenta with structure constants but such that the 
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corresponding Lie group is not compact. We will deal with this difficult case in a seperate paper 

In all the examples studied we recover the usual results. This might seem surprising because 
one would think, e.g. that the space of solutions to the single quadratic constraint M = p\ +P2 = 
is larger than the space of solutions to the individual linear constraints p\ = P2 = 0- Indeed, 
the general solution of the former is of the form fi(z) + f2(z) where z = x\ + 1x2 and /i,/2 
are smooth functions while the general solution to the latter are the constants where we have 
used the representation pj = ihd/dxj on 7i = L,2(R 2 ,d 2 x). However, solutions of the first type 
do not appear in the spectral resolution of the Master Constraint. Intuitively, this comes about 
because a physical Hilbert space based on square integrable linear combinations of holomorphic 
or antiholomorphic functions must be of the form L>z((D, dzdzp(\z\)) with a damping factor /o(|z|). 
However, this Hilbert space is not a representation space for a self adjoint representation of the 
Dirac obsrvables with respect to M: Indeed, the induced action of the Dirac observables pj is 
not self-adjoint in this representation. Thus, this representation is not viable unless we restrict 
ourselves to the constant functions. These are automatically selected by the spectral analysis 
of the Master Constraint which in turn is induced by the spectral analysis of the individual 
constraints Cj = pj . 

Another point worth mentioning is the following: Assume a simple situation, a finite number 
of first class constraints Cj which may close with structure functions only. In the case of non 
- trivial structure functions, the constraints Cj generically must not be self - adjoint for the 
following reason: Let {Cj, Ck} = fjk l Ci with non - trivial, real valued structure functions fjk 1 . 
Suppose that the Cj are self -adjoint then in quantum theory we expect a relation of the form 
[Cj,Ck] = ih(fjk l Ci + Ci(fjk where the symmetric ordering is forced on us due to the 

antisymmetry of the commutator. Since fjk 1 is real valued fjk 1 — (fjk should vanish as 
H — » 0. Now suppose that \lMs a generalized solution 6 of all constraints, Cj^f = for all j. 
Applying the commutator we find = ih[Ci, (fjk)]^ for all j, k. Now as h — > 0, the expression 
iH[Ci, (fjf.) ] becomes the Poisson bracket h 2 {Ci, fjk l } to lowest order in H. Unless this classical 
quantity vanishes, the solution \& not only satisfies Cj^> = but the additional constraints 
[Ci, (fjk) }^ = which have no classical counterpart. Iterating like this it could happen, in 
the worst case, that \E' satisfies an indefinite tower of new constraints, leaving us with the only 
solution \& = 0. In general the solution space will be too small as to capture the physics of 
the classical reduced phase space, displying a quantum anomaly. In the example with structure 
functions studied in this paper it was actually the case that {Q, fjk l } = which is why we 
succeeded in using self adjoint constraint operators there . This argument just given is not new, 
it was presented first in [21], but it is often forgotten. 

Finally, let us assume that zero lies in the point spectrum of all the Cj for a given Hilbert 
space representation and set M := YLjCjCj. Now the eigenvalue equations Cj^ = and 
M ^ = make sense for elements ^ £ 7i. Then we claim that that Cj^ = for all j if and 
only if M ^ = 0. The first implication is clear. For the converse we compute =< >= 
Y^j HC"^!! 2 - Thus we see that in this case the Master Constraint Operator captures exactly the 
same physics as the individual constraints, in particular it will suffer under potentially the same 
quantum anomalies. Hence, nothing is swept under the rug 8 . In case that zero also lies in the 
continuous spectrum, the DID procedure of [2] generalizes the argument just made. 



6 More precisely we should look for solutions in a space of distributions $* dual to a dense and invariant (under 
the constraints) subspace $ C H in the form ^((jjf) = for all / G $ and all j. 
7 We should study another example in the future where this is not the case. 

8 However, it is interesting to notice that by subtracting the minimum of the spectrum from the Master 
Constraint we may also be able to deal with anomalous situations! 
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